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Abstract 

The inflationary mechanism of mode amplification predicts that the state of each mode 
with a given wave vector is correlated to that of its partner mode with the opposite vector. 
This implies nonlocal correlations which leave their imprint on temperature anisotropies 
in the cosmic microwave background. Their spatial properties are best revealed by using 
local wave packets. This analysis shows that all density fluctuations giving rise the large 
scale structures originate in pairs which are born near the reheating. In fact each local 
density fluctuation is paired with an oppositely moving partner with opposite amplitude. 
To obtain these results we first apply a "wave packet transformation" with respect to one 
argument of the two point correlation function. A finer understanding of the correlations 
is then reached by making use of coherent states. The knowledge of the velocity field is 
required to extract the contribution of a single pair of wave packets. Otherwise, there is a 
two-folded degeneracy which gives three aligned wave packets arising from two pairs. The 
applicability of these methods to observational data is briefly discussed. 
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I. INTRODUCTION 

In inflationary models, primordial density fluctuations and primordial gravitational 
waves are described by Gaussian ensembles with well defined correlations in wave-vector 
space k. These correlations lead to the temporal coherence of the modes when re-entering 
the Hubble horizon in the adiabatic era [1-4]. As far as the primordial density fluctuations 
are concerned, the temporal coherence can be now considered as an observational fact since 
it is necessary to obtain multiple acoustic peaks in the spectrum of CMB anisotropics [4, 5]. 

So far the analysis of the correlations have been mostly performed in Fourier k-space, 
simply because they are diagonal in this representation. Nevertheless, it is of value to 
also analyse the two-point correlation function in the position representation. Indeed, 
as shown in [6] , this analysis displays the space-time causality of the mode amplification 
process. In this paper, we shall use a mixed representation based on local wave packets. 
This third analysis possesses its own virtues, and should be thought as complementary to 
the two other representations. In fact, we found it the most appropriate when focusing 
on the spatial correlations on a given time slice. This is because the use of local wave 
packets introduces spatial correlations by coupling different k-modes which were so far 
independent. Notice also that we shall work with the three dimensional Green function 
and not with its restriction to the Last Scattering Surface. This choice gives simpler 
expressions unencumbered by the projection on a 2-sphere. 

Wave packets are introduced by applying a "wave packet transformation" to the two- 
point correlation function with respect to one of its argument. Doing so one obtains 
a one-point function which displays a universal structure consisting of three local wave 
packets on a line. (To our knowledge this has not been noticed before). The central 
wave has relative amplitude two and corresponds to the chosen wave. The two others 
have amplitude minus one and correspond to two partners. As we shall later see, the 
reason for this three-folded structure follows from the fact that one deals with a snapshot 
of field configurations when using the two-point function on a given time slice. This 
implies that one cannot distinguish "left" moving from "right" moving configurations, 
thereby inducing a doubling of the partners. To raise the degeneracy, one needs to take 
into account the velocity of the waves, or equivalently distinguish positive from negative 
frequencies. Mathematically this raises no difficulty and can be obtained by using the 
Klein-Gordon product when applying the wave packet transform. Doing so, only two local 
wave packets of opposite amplitude are found, as one would have expected. 

There is a complementary way to interprete the use of wave packets: they can be viewed 
as introducing a filter in k-space. This opens the way to refine the procedure of filtering 
by working with the distribution of field configurations rather than with mean values. 
Indeed the Gaussian ensemble of field configurations determines, on one hand, the mean 
properties such as the two-point correlation function. These, together with Boltzmann 
equations [7, 8], determine the C\ the power of the temperature anisotropy multipoles. 
On the other hand, the knowledge of the ensemble also gives the probability to find a 
particular set of configurations, i.e. a particular realization of the ensemble. 

In the last two Sections, we exploit this second aspect in order to determine the spatial 
properties of the correlations associated with the realization of configurations described 
by a local wave packet. The procedure [9] consists in isolating these field configurations 
during the adiabatic era and to compute the correlations within that restricted set. These 
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correlations show up specific spatial properties which are some how smeared when dealing 
with the entire ensemble, i.e. with the mean values. These correlations have a double 
origin. First, their localization and their Fourier content depend on the chosen set. Sec- 
ond, their space-time structure is independent of this set and directly follows from the 
amplification process (or equivalently, the neglect of the decaying mode). 

It should be noticed that the space-time properties of these correlations coincide with 
those obtained by having applied a wave-packet transform to the two-point function. 
However, this second procedure is more general in that it gives also rise to correlations 
in amplitude. These cannot be obtained by working with the two-point function because 
in that case the mean has already been taken. Finally, even though the first procedure 
is simpler, the physical interpretation of the correlations it displays is unclear, at least to 
us. On the contrary, the interpretation of the analysis performed in configuration space 
is unambiguous and reached on a more fundamental level. The question of whether our 
procedures can be implemented to observational data is addressed at the end of the paper. 

II. BOGOLIUBOV TRANSFORMATION AND TWO-MODE STATES 

In this section we recall the basic elements which define Bogoliubov transformations 
in a the physical interpretation of the correlations it displays is unclear, at least to us. 
On the contrary, the interpretation of the analysis performed in configuration space is 
unambiguous and reached on a more fundamental level. The question of whether our 
procedures can be implemented to observational data is addressed at the mological context. 
In Eq. (9) we introduce the notion of two-mode states which will play a central role in 
encoding the correlations we shall focus on. 

It has been shown that the evolution of linearized cosmological perturbations (primor- 
dial gravitational waves and density perturbations) reduces to the propagation of real, 
massless, minimally coupled scalar fields in FRW spacetimes [3]. For simplicity, in this 
article, we shall study the fluctuating properties of a test scalar field £ in a homogeneous 
background. The translation of the results to physical fields represents no difficulty. 

We work with a line element with flat spacial surfaces: 



where Ti = d v a/a is the conformal Hubble parameter and V is the gradient with respect 
to the comoving coordinates x. It is convenient to introduce the rescaled field <p = a£ and 
to decompose it into Fourier modes 



ds 2 = a(r]) 2 [-drj 2 + 5ijdx l dx j ] . 



(1) 



The field £(77, x) obeys the d'Alembertian equation: 



<9^ + 2W^--V 2 £ = 0, 



(2) 




ikx 



(3) 



The time dependent mode 0k obeys 



{tf+u>l) k = O, 



(4) 
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where k = |k| and where the time-dependent frequency is given by 

^) = k 2 -^. (5) 
In second quantization, these modes are decomposed as 

4>k{v) = a^(pk(v) + aL k 4>* k (v) , (6) 

where the 'hat' characterizes operators. In this decomposition, 4>k is a solution of Eq. 
(4) with unit positive Wronskian [15]. It depends on the norm of k only since we work 
in an isotropic background. The operators and are the creation and annihilation 
operators of a quantum of comoving momentum k. The ground state of each mode k is 
defined by 

o k |0,k)=0. (7) 

Since annihilation operators of different momenta commute, the ground state of the field 
is the tensorial product over all k: 

|0> = n®|0,k>. (8) 

k 

When studying correlations due to pair creation, it is appropriate to re-write the vac- 
uum in terms of two-mode states: 

io) = n ® io, k> 2 . (9) 

k 

The tilde tensorial product takes into account only half of the modes and the k-th two- 
mode vacuum state is defined by 

|0, k) 2 = |0, k) ® |0, -k) . (10) 

(Notice that this notion can be generalized to other states whenever both modes are in 
the same 1-mode state.) When using this writing, one must pay attention not to count 
modes twice. To this end, for a real scalar field, one needs to separate (arbitrarily) the 
momentum space in two. For definiteness, we choose the separation according to the sign 
of k x , the x-component of the momentum. Then the product in Eq. (9) is performed 



1 The k-mode of the gravitational waves obeys Eq. (4) equation whereas the density fluctuation mode 
has a frequency given by 

2/ n 2, 2 Z 

z 

The function z(rj) is determined by the the background evolution: 

(H 2 -H') 1/2 
z = a- — '- . 

nc s 

During inflation c s = 1 and during the adiabatic era it is given by the sound velocity c 2 — Sp/Sp [3]. 
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over momenta with positive k x only. To emphasize this we shall call modes, states, and 
operators right (R) or left (L) according to the sign of k x . Hence we write = and 
d£ = a_k- Using this notation, the two-mode vacuum state obeys 

ag\0, k) 2 = a£|0, k) 2 = 0. (11) 

In non-stationary backgrounds, the frequency Eq. (5) depends on time. Hence the non- 
adiabaticity of the propagation leads to spontaneous excitations of the various modes. To 
characterize these transitions, it is appropriate to introduce two sets of modes. These are 
positive frequency solutions of Eq. (4) at early and late time. In Appendix A, they are 
explicitly given when considering a cosmological evolution which starts with an inflation- 
ary phase and ends by a matter dominated period after having experienced a radiation 
dominated period. As usual, we shall use the labels 'in' and 'out' to designate states and 
operators which are defined with respect to the corresponding modes. Since Eq. (4) is 
homogeneous and linear, in and out modes are related by a Bogoliubov transformation 

cw = «^r t w+^c t, ('?)- (12) 

The corresponding transformation between in and out operators is 

~R,in * ~R,out n ~L,outi /io\ 

a k = «fc a k "A a k ■ ( 13 ) 

Because of the homogeneity of the background, this transformation is 2 x 2 block-diagonal 
as it couples k to — k only. Hence every produced k-particle will be accompanied by a 
partner of momentum — k. Moreover particles characterized by different momenta are 
incoherent in the in vacuum (in the sense that in the expectation value of any product of 
annihilation and creation operators of different momenta will factorize). 

These two properties are made explicit when expressing the in vacuum in terms of out 
states (i.e. states with a definite out particle content). From Eq. (13), using the notations 
of Eq. (9), one gets (see App. B in [9], [16]) 

|0, in) = Y{ ® |0, k,m) 2 

k 

= n^(^ ex p(^«k*° Utt «k OU * t )|0,k,ont) 2 ) , (14) 

where = (3k/a* k . From this writing we see that the in vacuum factorizes into a product 
over half the momenta of sums of two-mode out states. It has to be emphasized that these 
out states carry no 3-momentum since they contain exactly the same number of R and L 
out k-particles. 

Our aim is to analyze how the these properties determine the space-time structure of 
the correlations of the <j> field. We shall use two different approaches. In Section III and 
IV we shall work directly with the two-point correlation function 

G in (r],x.;r)',x!) = (Om|^(r/,x) <£(t/, x')|0m) . (15) 

In Section V and VI, we develop an alternative approach, based on [9], which is more 
fundamental as it is based on the correlations in configuration space encoded in Eq. (14). 
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III. SPATIAL CORRELATIONS INDUCED BY PAIR-CREATION 

Since we are dealing with a free field, all (in-in) expectation values of products of the 
field operator can be decomposed in terms of the two-point function Gi n . Its late time 
properties are best revealed by decomposing the field operator into out modes. One gets 

Gi n (r),x;r]' ',x') = G out (r], x; rj' , x') 

+2 l^k^^nkRe^r^r*^') + ^<T^)CV)} • (16) 

In the first line, G ou t is the Wightman function evaluated in the out vacuum. This quantum 
contribution is 0(1) whereas the second term is proportional to the occupation number 
n k . Hence when n k 3> 1 the vacuum contribution can be neglected (unless one computes 
operators containing commutators, because in certain cases, the second term might not 
contribute since it is symmetric in x,x'). Notice also that we could have split Gi n into a 
commutator and an anti-commutator. In the large occupation number limit, the dominant 
terms coincide. 

The second line of Eq. (16) is governed by two quantities. First one has a diagonal 
term 

Kr K'tRhn = Kl «^L>in = U k 5 3 (k - k') = \a k \ 2 \z k \ 2 5 3 (k - k') , (17) 

which fixes the mean number n k . (The symbol ( • ); n designates the in vacuum expectation 
value: (Om| • |0m).) Second one has an interfering term 

«fl «k" i)in = c k S 3 (k - k') = \a k \ 2 z k 5 3 (k - k') , (18) 

which governs the coherence of the distribution. By coherence we mean that the expec- 
tation value of a product does not factorize. In the present case, it is ^k"!)" 1 ^ 
in in = which thus expresses the coherence. For incoherent distributions, 
such as thermal baths, one would get c k = for all k. Notice also that the expec- 
tation values which differ from the above ones by one additional f on an operator a 
all vanish. This last property is valid for all homogeneous and isotropic distributions 
(and not only those resulting from pair creation). The degree of coherence of the dis- 
tribution is given 2 by \c k \/(n k + 1/2) < 1. For pair creation from vacuum, one has 
\c k \/(n k + 1/2) = 2\z k \/(l + \z k \ 2 ). 

For macroscopic occupation numbers, the norms of the diagonal and interfering terms 
coincide since \z k \ = \(3 k /a k \ — > 1. In this limit the interfering term can thus be written as 

c k = -n k e i2 ^. (19) 

Taking into account the isotropy of the distribution, the dominant part of the two-point 
function simplifies and reads 

G»fex;i',x') = f ^ ! l^p ^ /ra { e *^W} /™ {e*^^)} (20) 



2 In this definition, n*, + 1/2 and Ck are the non- vanishing elements of the covariance matrix in the two- 
mode state |0, k, 171)2, defined by the expectation values of the anticommutators of clr, &l, and a\. 
This matrix is also the covariance matrix of the corresponding classical distribution [14] . 
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The integrand is a product of two classical waves (= Im { e l ^ k (f> k mt } ) . Three remarks are 
in order. First, this factorization could not have been performed if the distribution did not 
obey \ck\ = n k . Therefore the fact that it can be done is an expression of the two- mode 
coherence of the underlying distribution. 

Second, when working with a (two-mode) coherent state, the Green function is also a 
product of two classical waves, see Eq. (C4). This therefore suggests to view Eq. (20) as 
resulting from an ensemble of two-mode coherent states, see Section V. 

Third, the two-mode coherence giving rise to the classical waves im{e*0£ u *} yields 
the usual description [10-13] based on the neglect of the decaying mode, see [14] for more 
details. For simplicity, let us consider a radiation dominated universe which follows a 
period of inflation. In this case the classical wave with kr/ r <C 1, where r/ r is the time of 
reheating defined in Eq. (Al), are proportional to s'm(kAr]) up to a correction term of 
the order of (fcr/ r ) 3 ~ 10~ 75 when inflation lasts for about 60 e-folds. They correspond to 
growing modes since the conformal time-lapse At? is proportional to a(rf). Given this strict 
correspondence between Im { e l ^ k } and growing modes, one can abandon the quantum 
settings and proceed with the effective description based on growing modes with stochastic 
amplitudes. In this paper we shall nevertheless use the quantum formalism in Sections V 
and VI for the following reason. It allows to treat separately right and left moving wave 
packets, a feature which is useful and which leads to a transparent interpretation of the 
results. The next Section however is based on the two-point function and can therefore 
be interpreted in either formalism. 

We now wish to illustrate how the above coherence (or equivalently the neglect of the 
decaying mode) induces spatial structures on a given time slice, e.g. on the Last Scattering 
Surface (LSS). The lapse Ar] appearing in the classical waves determines the characteristic 
size of the structures on the LSS. Indeed when rj' = rj, the stationary phase condition 
applied to the integrand of Eq. (20) gives two solutions. First one gets |x — x'| = which 
is responsible for the usual divergence in coincidence point limit. More importantly, there 
also exists a non-trivial solution: 

|x - x'| = 2At] = 2( V - 8 k 4> k ) = 2( V - 2r] r ) . (21) 

This only results from the interference term c&. The lapse 2Ar] designates the mean 
separation reached by the particles and their partners from their birth near a = 0, see 
Appendix A. This interpretation will become clear in the sequel. 

It should be also noticed that the contribution of this second term is negative, thereby 
causing a dip in the two-point function [6]. The origin of this dip can be traced to Eq. 
(19) and the fact that V'fe 1. It tells us that only the growing mode has been kept. As 
we shall demonstrate in Section VI, this implies that the partner of any local over-density 
is a local under-density. Notice finally that the relative weight of the usual solution and 
that of Eq. (21) is two. In the next Section this factor shall be recovered and explained 
in terms of local wave packets. 

IV. WAVE PACKET TRANSFORM AND SPACE-TIME CORRELATIONS 

In this section we analyze the space-time correlations obtained by wave packet trans- 
forming the two-point function with respect to a i?-moving wave packet. We shall consider 
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two different scalar products: one based on the usual overlap and the other based on the 
Klein-Gordon product. The first product leads to a spatial structures containing three 
packets whereas the second leads to two packets only. The reason comes from the fact 
that .R-moving and L-moving modes equally contribute in the first case, thereby leading 
to a doubling of the partners. The knowledge of the velocity field is required to lift this 
degeneracy. 

A i?-moving wave packet can be written as 

I Til I ' 

OR, " 



, Vfl (r?,x) = jd*k ^J^#*fo) + C . c .) 



= 4tk + 4:1, • (22) 

where the tilde integral means that only positive values of k x are considered. The symbol 
Vr designates the set of Fourier amplitudes v k- It is used to remind that the specification of 
the wave packet has been made in the i?-sector. In the second line, we have decomposed 
the wave into its positive and negative frequency content. This will be needed when 
considering the Klein-Gordon product. 

To be specific we shall consider a single Gaussian wave packet. For clarity of the 
equations, we write 

Vk = vf R (k), (23) 
where v is real and positive and where the function is normalized to unity: 

'<Pk\f R \ 2 = l. (24) 

We shall use the following Gaussian wave 



f R (k) = e^Ne~ I ^e- ikx V fo?0 , (25) 

where N > is a constant such that Eq. (24) is satisfied. The mean momentum of the 
wave packet is k and its mean position at r]o is xo. 4> is the phase of the positive frequency 
part of the wave evaluated at xq,t]o. (Notice that when considering only the 1-particle 
sector, this phase would be inaccessible. Instead when dealing with coherent states, it 
is an observable). Suppose also that we are working in radiation dominated era, so that 
positive frequency modes are 



'k 



1 

2k 



? M *(?7) = -= e - <fe » . (26) 



The i?-wave is thus 



4) RVr = vN [ d?k-^e- {± ^- ( e *(x-xo) e -ifc(r;-»») e ^ + C A . (27) 
J V2A; ^ ' 

By making use of the stationary phase condition, one finds, as expected, that this wave- 
packet is maximum along the classical light-like trajectory 

XR(r/,k) = x + (77 - r? )lk, (28) 

which passes through xo at t]q with a momentum k. (The vector l k is the unit vector in 
the direction of the velocity of the chosen wave.) 
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A. Wave packet transform based on the usual product 

We use the above i?-moving wave to "wave-packet transform" the two-point function 
with respect to one of its arguments: 

0(x,t/;Vr) = Jd 3 x' 4> R (i] ,x)Gi n (ri,x;ri ,x) . (29) 

In Fourier transform we get 

e(x,r ? ;V ii ) = Jj^L. (v k ^r(Vo)+cc) CWC*(%). (30) 

By keeping only the growing mode, see Eq. (A7), and using Eq. (26), one obtains 

e(x,r ? ;Vi?) = ^ Vr (x-xo,i]-i] ) + ^ Vr (x-x ,-(i]-i]o)) 

-^V fl (x - x , -(r? + % - 47/y)) - ^v fl (x - X ,7/ + 7/0 - 4?7 r ) . (31) 

The function ^v fl is defined by 
*v«(x,r/) = 

= * « + $H j (32) 

where we have decomposed the amplitude of Eq. (25) as its norm times its phase in 
order to exhibit the linear phases in x and rj. In Eq. (31) it is the same function ^v R which 
governs the four contributions. The reason for this simplicity arises from the fact that 
the differences between the integrands in Eq. (30) are given by Zk = flk/oi* k whose phase 
is linear in k, see Eq. (A6). Hence the four contributions only differ by their temporal 
argument and their relative sign. This relative sign guarantees that the integral over x of 
vanishes irrespectively of the chosen wave packet. 

It is of interest to analyze the case where the integral Eq. (32) can be evaluated by a 
saddle-point approximation. Then the shifts in r\ translates into shifts in x. Indeed, given 
a wave-packet of mean momentum k, Eq. (31) becomes 



e(x,7/;V R ) = $ Vfl (x-x i? (7/,k); E ±tR ) + $ Vfl (x - x^t/, -k); Ej. )fl ) 

-$v R (x-x L (7/,k); £ ±jL ) - $y H (x-x L (7/,-k); E ±jL ) . 



(33) 



The spreads £j_ 5 _r and ^ are rather complicate matrices (in 3D). We postpone the 
analysis of their interesting properties in a separate subsection. The function xr has 
already been defined in Eq. (28). Notice that we have enlarged its definition to "negative" 
wave vector so as to describe the second term of Eq. (33) which is a L- moving wave packet. 
The new function xl is defined by 

XL(7/,k) = X i? (-7/ + 47/ r ,k), 

= x - (7/ + 7/o-47/ r )l S . (34) 
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It corresponds to the trajectory of the partner of the wave of momentum k. This will be 
clearly established in Section VI when dealing with coherent states. We can already notice 
that the separation between the two waves is 

x r (t}, k)-x L (r?,k) = 2(77-27?,) l t . (35) 

The interpretation of this result is clear: it is the separation reached by the two waves 
since their creation (amplification) near the big bang. Indeed, rj — 2r/ r oc a{rj) during the 
radiation dominated era, see Appendix A. 

The distance \xr — x/J is universal in the following sense. First, as expected, it is 
independent of both xo (because of the homogeneity of the process) and the direction 
specified by 1^ (because of isotropy). Somehow more surprisingly 3 it is also independent 
of the norm of k. This results from the conformal (or scale) invariance of the theory. This 
independence of the traveled distance is therefore complementary to the well known fact 
that the power spectrum is (nearly) scale invariant. When considering the field at the 
origin of primordial density fluctuations, see e.g. [25], this result implies that all wave 
packets at the origin of the large scale structures are born near the reheating time, see 
Figure 3. The same remark also applies to primordial gravity waves. 

It should also be remarked that Eq. (35) is the version of Eq. (21) wherein one has fixed 
the direction specified by the mean momentum k. The role of the wave packet transform 
is therefore to isolate from the mean the contribution specified by 4>r, i.e. by the set of 
Fourier components t>k- 

We now proceed with the description of the space-time properties of four waves in 
G(x, 7/; Vr). To this end we first represent in Figure 1 the function for a one dimensional 
wave packet. By construction the two wave packets governed by x#(±k) merge at rj = tjq, 
see Eq. (31) and Eq. (33). Thus, at that time we have a three-folded picture: a central 
wave with weight 2 (since both right and left movers equally contribute) surrounded by 
the two partner waves of weight 1 which are respectively a L and a R moving waves. These 
are shifted to the left and to the right with respect to the central wave in the direction 
specified by k. Their maxima are located along the classical trajectories Eq. (34). One 
clearly sees the factor 2 in amplitude and the phase opposition of the partner waves. One 
also verifies that the integral J dx@(rj, x) identically vanishes for all wave packets. 

We thus see that the wave packet transform with respect to a i?-moving wave has 
isolated two pairs of wave packets. The doubling arises from the fact that the Fourier 
transform in Eq. (30) is insensitive to the velocity of the wave. Hence the wave packets 
with mean momentum ±k equally contribute to 0(7/, x). This two-folded degeneracy can 
be lifted if one works with a product which is sensitive to the velocity of the waves. To 
this end we shall use the Klein-Gordon product in the next subsection. 



3 This is at least not usual in that this is not what is found when considering pair production of massive 
particles in de Sitter space, eletro-production in a constant electric field, or pair production giving rise 
to Hawking flux in quantum black hole physics. In all those cases, the norm of the wave vector does 
characterize pairs produced at different times, see [9, 16]. 
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FIG. 1: On the left, we have plotted 6(77, x) of Eq. (29) in a 1 dimensional case. We have rescalled 
the spatial direction by c s the speed of the wave. The wave-packets have a mean momentum 
c s k — 7r/3, a spread a = 3 and a phase = 0. The chosen R-wave is centered on x = at 770 = 3. 
Reheating is at i] r = — 1CP 2 . (Hence one has c s kr]o = w, i.e. k corresponds to the first peak in 
the CMB if 770 corresponds to the recombination.) The conformal time ranges from reheating to 
i] = r)o . On the right, we present a spatial section at 77 = 7) . 



B. The Klein-Gordon product 

When using the Klein-Gordon scalar product in the place of the simple product of Eq. 
(29), one gets 

e^ G (x,r?; V R ) = |dVG m (r),x,))o,x' )id m - c.c.) = 

/: pikx 

= 2 J Re(^^+ ) (x-x , 77-770) 

-tfytfg^x - x , -(r) + r/o - 4r/ r )) , (36) 

where the positive frequency wave is defined in Eq. (32). On the first line, we have 
used the difference between the positive and negative frequency components of <pR in order 
to cancel the minus sign which appears in the Klein Gordon product. On the second line 
we thus get the real part of the integrand. It differs from that of Eqs. (31,32) in two 
respects. First there is an extra factor of k which arises from the derivative with respect 
to 77. Second and more importantly, the KG product of wave with opposite frequencies (i.e. 
velocities for a given k) vanishes. Hence this leads to a reduction of the four contributions 
of Eq. (31) to two waves forming a single pair. 

When using the KG product between the two-point function Gi n and a i?-moving wave, 
one thus correctly isolates the pair whose i?-moving mode corresponds to the chosen wave. 
Together with Eq. (31), this is the main result of this Section. 
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The extraction of the contribution of a single pair is clearly displayed in Figure 2. 
The two wave packets originate near the reheating from a small patch which is centered 
at (xfl + x/,)/2. In fact the two waves travel on the would be particle horizon, i.e. the 
particle horizon if there were no inflation before rj r . This can be seen from the fact that 
the solution of — x^ = 0, i.e. the tip of the lightcone, gives rj = —2rj r which would 
correspond to a = 0. The spatial extention of the patch at the reheating is governed by the 
spread of the waves: 1/a. From these results, one finds that in a radiation dominated era, 
the partner of a gravitational wave is always outside the Hubble radius centered on the 
other wave. This means that the coherence will never be detectable by any measurement 
performed within a Hubble radius. (In this one gets a situation similar to that of Hawking 
radiation since the partners of Hawking quanta are all inside the horizon[16]). However, 
thanks to the imprint they leave on the LSS, both members can now be seen by us (when 
they are properly aligned, see Figure 3). 

As previously discussed, the two waves are in phase opposition. Moreover, because of 
scale invariance, the above properties are valid for all wave packets. 




FIG. 2: On the left, wc have plotted &kg in conformal coordinates. The conventions are identical 
to those of Figure 1. On the right we have presented sections at r\ = 770 = 3 of this amplitude 
for two different choices of the phase </>: <fi = 0, tt/3. In each case, the partner wave is in phase 
opposition in the sense that Eq. (59). The two waves are separated by 2 (770 — 2r] r ) ~ 6. This 
filtering giving rise to the contribution of a single pair will be recovered in Section VI by performing 
a projection in configuration space rather than a product in space-time as in Eq. (36). 



C. Spreads 

It is of great interest to analyze the spreads of the various contributions in Eq. (33) or 
Eq. (36). Indeed their properties can be exploited to reveal the presence of the correlations 
between R and L sectors. 
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The result of the saddle-point evaluation of Eq. (32) is given by, see [9]: 



$ Vfl (x) ~ Re 



=e x(k) exp ^^.(s-^.x/) 



(37) 



where the function x(k) regroups the terms in the exponential of the integrand in Eq. (32). 
The matrix £ is defined by 2£jj = — didjxlk where i,j = 1, 2, 3 label the three conformal 
coordinates. The three eigenvalues of £jj give the spreads in position of the Gaussian 
wave-packet Eq. (37). In the case of a single wave-packet with axial symmetry around 
k, the three eigenvalues reduce to two scalars. The first one is the spread in direction 
longitudinal to k. It is independent of time and equal to 1/|<t|, as in one dimension. The 
other governs the spread in the directions orthogonal to k. It grows linearly with time, 
as in the case of a non-relativistic wave-packet. Its value therefore differs for the R and L 
wave: 

2^(7?) = ^+^°, (38a) 

2S ± , L (r?) = ^2 + l 1 • ( 38b ) 

When a is taken real, the spread of the R-wave is minimal at rj = 770, as expected. 
(Remember that the spread in the perpendicular directions is given by the square root of 
the norm of £). The spread of the partner L-wave is thus larger as it is governed by the 
accumulated conformal time from the detection of the R wave at t/o, back to rj r where the 
pair emerges, and then forward to 770, for more details see [9]. The relative increase of 
L-wave spread becomes large for short wave lengths modes. As an example, for a wave- 
packet of mean momentum k(rjo — 2r/ r ) = 10 (this would correspond to the third peak in 
the CMB anisotropies spectrum) and spread a = k, one has |£/_,(77o)/£.r(77o)| ~ 6. 

It is now crucial to notice that one can fine tune the imaginary part of a to obtain a 
partner wave which is more peaked than the wave which is chosen in Eq. (29). Indeed, 
one can exploit the coherence of the wave packet and that of the vacuum, i.e. the k 
dependence of Eqs. (40, A6), to obtain constructive interferences around by canceling 
the accumulated effect due to the total time of flight. The maximal effect, i.e. the minimal 
norm of £1,(7/0), is reached when taking Im(o-~ 2 ) = — 4(7/0 — 2rj r )/k. When considering 
the above example with k(rjo — 2rj r ) = 10, one gets 

^(t/oV^MI-^, (39) 

that is, the perpendicular spread in space of the partner wave is about six times smaller 
than that of the wave which is used in the integration of Eq. (36). This reduction of the 
partner wave spread is a direct consequence of the two-mode R — L coherence. Hence it 
can be used as a check to quantify the degree of coherence given some observational data. 
It should be also noticed that it is not necessary to identify the velocity in order to bring 
this effect into evidence. Indeed this reduction equally applies to the spreads of the last 
two waves in Eq. (31). (In that case, the above symbols R and L should be interpreted 
in their generalized sense given after Eq. (33).) 

In the next Sections, we present a more sophisticated procedure to extract space- 
time correlations. It is based on the introduction of a filter in configuration space before 
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computing expectation values. Being not based on mean values, it allows to generalyse 
the former analysis. In particular it reveals correlations in the amplitude of wave packets. 
The reader not interested by these developments can proceed directly to the Conclusions 
where he will find a resume of the results. 



V. COHERENT STATES, CLASSICAL WAVES AND RANDOM PROCESSES 

We introduce coherent states of the field 4>. They are the key elements for our analysis of 
the correlations associated with a realization of a particular set of classical configurations. 
In quantum terms such a realization can be described by the detection of the corresponding 
states, namely coherent states. Moreover, since we want to select localized waves, we need 
to form wave packets by summing different k-modes. Therefore, we shall work with 
coherent states of these wave packets. 4 

The filtering of a set of classical waves is implemented by introducing a projector on the 
corresponding coherent state. The introduction of this projector modifies the correlations 
which existed in the "in vacuum", i.e. the state of the system which encodes the pair 
creation process. New correlations are introduced at the expense of reducing the pre- 
existing ones. In particular spatial correlations are generated by detecting a local wave. 



A. EPR correlations and conditional values 

To describe the correlations which exist in the in vacuum, we analyze correlations 
amongst out states. To justify this analysis in the present context, let's consider the 
following gedanken experiment. Suppose that n i?-out-particles of momentum k have 
been detected and that nothing is known for the other modes except that, before this 
detection, the field was in the in vacuum. Given that detection, one can ask what are the 
probabilities to find m particles of momentum k'. The answers to this type of questions 
are governed by the partially 'reduced' state obtained by projecting the in vacuum onto 
the state characterizing the (partial) detection, namely |n, k, R, out): 

( Zu) n ■ 

l r ed n k _r) = (n, k, R, out\0, 111)2 = — r \n, k, L, out) \\ <£> |0, k', in) 2 ■ (40) 

\ctb\ ■ LJ - 

In the k-th two-mode sector, one finds the one-mode L-state entangled to \n, h,R,out). 
It is a pure state with the same occupation number. This results from the EPR-type 



4 The use of coherent states can be conceived from several point of views: either as a mathematical way 
to introduce classical waves in quantum terms, or more physically, as resulting from a detection of such 
waves, or even more intrinsically, through decoherence induced by interactions. Indeed oscillators weakly 
coupled to an environment evolve into coherent states [17, 18]. Therefore it is to be expected that in 
cosmology the weak non-linearities which are generally ignored will replace the pure (two-mode squeezed) 
state by a mixture of (two-mode) coherent states [14]. However, to our knowledge, the properties of 
this decoherence process (when it occurs, how it modifies the state of the k-modes, and how to put it 
into evidence) have not been fully derived. Finally, the use of coherent states provides an interesting 
alternative to [10-13] when analysing the emergence of classical and stochastic properties in inflation. 
Indeed the detection of coherent states is a Gaussian random process. Notice also that coherent states 
have been already used in [24, 25] to study the semi-classical limit. 
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correlations present in the in vacuum, see Eq. (14). (These correlations deserve the 
label 'EPR' since they are of the same character as those encountered with spins: the 
spin projection on a axis is here replaced by the occupation number. Notice also that in 
both cases, a symmetry is at the origin of the entanglement, rotation invariance there, 
translation invariance here.) The two-mode sectors with k' 7^ k are all unaffected by the 
detection. Therefore the probabilities to find particles of momentum k' / k are unchanged. 

More generally, according to the rules of quantum mechanics, the "conditional" expec- 
tation values which result from the above detection can be written as 

(O)-^, (41) 

\ll/in 

where the projector onto the detected state |n, h,R,out) is 

tl n ,k,R,out = \n,k, R, out)(n, k, R, out\ O l kiL ® lkVk • (42) 

In Eq. (41) we have used ft 2 = ft to simplify the denominator. 5 The above projector 
is "partial" in that it is unity in all sectors but the iZ-mode k. 

The use of projectors which act only on the R sector is particularly interesting because it 
guarantees that the L-mode content will be only determined by the correlations which are 
present in the state of the field. It is therefore the appropriate tool to unravel the intrinsic 
properties of the R — L correlations. The next step is to determine how these correlations 
give rise to the spatial properties of conditional values. This is done in Sections VI. Notice 
finally that we shall no longer use projectors which specify the occupation number as in 
Eq. (42). Instead we shall use projectors based on coherent states because these properly 
characterize the various realizations of the ensemble, see footnote 3. 



B. Coherent states and classical waves 

Because of the entanglement between R and L modes in the in vacuum, only half of the 
modes are independent. Hence we shall use coherent states which are defined in the right 
sector only. We also work with out modes because the detection of waves is performed at 
late time, i.e. during the radiation or the matter dominated epoch. In this Section, we 
consider coherent states which specify the amplitudes of all iZ-modes: 

\V,R) =]\®\v k ,R), (43) 

k 



5 A technical comment is in order for the readers familiar with the work of Aharonov et al. [19] or with its 
applications to pair creation [9, 16, 20]. In these works, a different conditional value of O, called weak 
value, was used. It is given by (II£>) / (It) in the place of Eq. (41). Mathematically, the difference between 
the two expressions is that weak values are generically complex whereas Eq. (41) is real for hermitian 
operators. On the other hand, when O and II commute they coincide since II 2 = IT. To understand 
the physical relevance of these two versions is a subtle question: which version should one use in a 
given context when some information concerning the final outcome is known ? The reader interested by 
this type of question will consult the original references. In inflationary cosmology, because of the high 
occupation number, a simplification occurs: the differences between the two versions are subdominant 
(i.e. 0(l/n)) since they arise from commutators. 



16 



where l^k, -R) is a (one-mode) coherent state of complex amplitude v^, see Appendix B. 

The product \V,R) is an eigenstate of 4>^\ the right-moving positive frequency part 
of the field operator : 



LOUtl 



= &W,R). (44) 

where the tilde on d 3 k means that one integrates over R modes only, i.e. k x is integrated 
from to oo. The function cf>^ is complex because it contains only positive frequencies. 
To get a real function, one should consider the expectation value of the observable 4>r, the 
field amplitude itself: 

4> R = (V,R\]> R \V,R) = (V,R0 R +) + ^ R ] )\V,R) = 2Re[^ ] ] . (45) 

The same can be done to the conjugate momentum of the field, see Eq. (B6). 

We can now verify that the (dominant part of the) expectation values computed in the 
state | V, R) can be expressed directly in terms of the mean value 4>r, in the same way that 
Eq. (20) was expressed in terms of the two sine functions. 

Using the normalisation of Eq. (24), the amplitude v is related to the expectation value 
of the occupation number 

(V,R\N R \V,R) = JcPk(V,R\a^ag\V,R) = v 2 . (46) 



It is also related to the current of </>^ : 



d 6 x <4 _) id v 4> ( r ] = v 2 . (47) 

Thus one has two alternative descriptions. Either one uses the quantum description in 
terms of the expectation value based on the counting operator or the classical concept of 
current based on the mean wave ■ The same conclusion is valid for other quantities 
such as Green functions, see Appendix B, or the 3-momentum. In all cases, when v 2 S> 1, 
the classical expressions based of the mean field cpR coincide with the corresponding ex- 
pectation values evaluated in the coherent state \V,R). The reason is that the ambiguities 
of operator ordering lead to differences governed by commutators which are subdominant 
in the large occupation number limit. 



C. Detections and random processes 

In this subsection we show two important results. First, when the (Heisenberg) state is 
the in-vacuum, the detection of the i?-moving configurations 4>R described by the coherent 
state |V, R) is a stochastic Gaussian process. This result is exact: it requires no approx- 
imation and is valid even before applying any decoherence process. Notice however that, 
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because of the entanglement between the R and L sectors, the particle content of only- 
half the modes (e.g. the i?-moving configurations) should be specified to get this result. 
Second, the detection of (pR fixes the L modes to be also described by coherent states. 
This result is also exact and follows from the Gaussianity of squeezed states and the R — L 
entanglement. These two results entirely determine the correlations which result from the 
detection of semi-classical configurations described by coherent states. 

To determine the consequences of a detection, it is appropriate to introduce the asso- 
ciated projector, see subsection VA. In the present case, it is 

U Vr = \V,R)(V,R\®1 l . (48) 

It is non trivial in the .R-sector only. The probability to detect the classical wave 4>r is 
given by 

P^ R = (0m|n Vfl |0m) = nK n | 2 , (49) 

k 

where the amplitude for the k-mode is 

h^L 2 



|«fc| 



, (50) 



see Eq. (C7) for the details. The probability Eq. (49) defines a normalized gaussian 
distribution for each i?-mode in every two-mode sectors. The normalization follows from 
the "density" of coherent states, see Eq. (C8). As already mentioned, only i?-modes have 
been so far specified. Had we performed a projection on both right and left sectors, the 
probability would have been exponentially smaller. Indeed, the ratio of the probabilities 
with or without double projection is, see Eqs. (52) and (Bll), 

pin 

^^=n-p(-K-^ k i 2 )' (si) 

r v R k 

where is the amplitude of the coherent wave in the L-sector. 

This exponentially suppression results from the entanglement between the left and right 
sectors. Indeed, when applying the projector fly^ on the in vacuum one gets 

n Vfl |0,m) = Y[A™\v k ,k,R)®\z k vl,k,L) . (52) 

k 

It is remarkable that in each two-mode sector, the L-state is also a coherent state. This 
results from Eq. (40), see also Eq. (C6). The L-mode amplitude is Zfct>£. It is fixed by the 
R- amplitude and by the pair creation process which is governed by Zk- The properties 
of the space time patterns we shall later exhibit directly follow from this double origin. 

Taken together, Eqs. (49-52) show that the notion of stochastic processes naturally 
emerge when questioning the in vacuum by making use of coherent out states. More 
precisely we have the following. Firstly, as one might have expected, the i?-mode amplitude 
i>k is a Gaussian stochastic variable of variance equal to \a^\ = n& + 1 ~ n^. Secondly, 
the k-th L-mode amplitude is "slave driven" by the detection of the i?-mode in that its 
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probability is centered around z k v k with a spread equal to 1, see Eq. (51). Therefore, 
in the large n k limit, this spread is negligible and one can consider that the L-mode 
amplitude is equal to z k v k . Thus, in the stochastic description as well, P v ™ is a two- mode 
distribution, as clearly seen from Eq. (52). 

These properties offer an alternative way to express expectation values in the in vacuum. 
It suffices to apply the following substitution: first, at the level of amplitudes d^ — > t>k, 
d£ — > z k v k , and second at the level of the distributions, the quantum distribution pi n = 
|0m)(0m| should be replaced by Py™ of Eq. (49). Notice that no dynamical assumption 
was needed, nor was it necessary to follow the time evolution of the modes. 

The emergence of classicity rests on the high occupation number n k S> 1 and on the 
(restricted) set of questions formed by inquiring about the coherent out state content of 
the in vacuum. Even though these conclusions are not new [10-12], the derivation which 
makes use of coherent states is particularly clear. In particular, it disentangles the question 
of the late time description of the in vacuum in the above stochastic terms from the more 
difficult question which concerns the evaluation of the time from which this stochastic 
description is valid. This time is determined by the efficiency of decoherence processes in 
the early cosmology, a subject not addressed in the present paper [13, 18]. 



D. 1-point and 2-point functions 

Besides the above substitution, one can consider the projector fly fl as in subsection 
V A, namely as defining a new ensemble of configurations with modified expectation values 
given by Eq. (41). It is of interest to present these expectation values in some details. 
Starting with 1-point functions, we have 

<^fl>V« = «k, (&Z?L)v R = Zk<- (53) 

In the in vacuum we had (a^'fym = (^"Dm = 0- The interpretation of the modification is 
clear: once we know that the classical wave 4>r has been detected, the mean R- amplitudes 
of the R k-modes are those of that wave. Moreover because of the EPR correlations in 
the in vacuum, the mean amplitudes of the associated L-modes are fixed by the detection 
of the i?-wave <j)R and z k . 

For the 2-point functions we have, 

<<°S <,r)v r = i*V , (4T K'!l)v r = z* k v k z k ,v* k , , (54a) 
(ai: t R a^ L )v R =v k z k/ v* kl , (54b) 
("i™R* K,L)v R = vUk<>. (54c) 

In the first line, the main modification with respect to in vacuum correlations is the loss 
of the diagonal character in k. This radical change follows from the strength of the pro- 
jection induced by II y n . Since all i?-k components are now described by coherent states, 
the above 2-point functions are entirely given by a disconnected contribution. For these 
2-point functions, the correspondence mentioned in section IV. B is exact. However, this 
is not the case in general because of non-vanishing commutators (consider for instance 
(^k"i? Q'^rIVr)- It * s onr y i n the large occupation number regime that the operator or- 
dering gives subdominant corrections. In the second line, we see that the in vacuum 
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correlations between R and L modes have been replaced by the "coherent state correla- 
tions" described by Eq. (52). They fix (a^ ut L ) VR in terms of {a^ R )v R . Notice finally that 
the variances of ano - ^k"*L vamsn > see Eq. (B4). The detection of a coherent state of 
the field can thus be seen as providing one classical realization of the stochastic ensemble. 

Before examining the correlations in space induced by the detection of 4>r, it is of value 
to determine to what extend one recovers (in vacuum) mean values from these conditional 
expectation values and from the distribution Py™ . Using Eq. (51) and Eq. (50), one finds 
that the ensemble average is defined by 



JVv p V W(1 P^ Wl 

z fc K| 2 <5 3 (k-k'). (55) 



JVv p v^Zk'vl, P™ R _ 



Pv p P v « 

The tilde over the functional integration is there to remind that the integration variables v k 
(%) are defined only for k x > 0, thus Vv^ = n^k- The result of Eq. (55) is in agreement 

k 

with the quantum result Eq. (18). For the diagonal R term, the correspondence between 
the ensemble average and the quantum result Eq. (17) is not exact and the mean values 
differ by a factor equal to |^| 2 . The origin of this discrepancy is that the R ctk',.R 

does not commute with the projector tly R , see also the third footnote. However, since 
the discrepancy rests on commutators, the two versions will agree when — ► oo. This 
agreement in the large occupation number limit confirms that field configurations are 
effectively characterized by a set of stochastic variables V with the two-mode probability 
distribution P in (V). 



VI. SPATIAL CORRELATIONS 

In the preceding subsection we gave the new expectation values when having detected 
the classical right moving configuration (f>R. Here we shall see that this detection leads 
to specific correlations in space. To exhibit these correlations it suffices to compute the 
modified expectation value of the field amplitude: 

(4>{v, x)) Vjl = 4>R, v R + <Pl, v R , (56) 
where the R and L mean waves are 

4>R,V R {r},*) = j&k(v\i J^^^k Ut (v) + cc)j , (57a) 

h,V R (r),x) = JM (z k vt ±^^\ri) + c.c.) . (57b) 

Eq. (57a) should cause no surprise. Since the projector tiv R completely specifies R- 
configurations, expectation values in the R sector are given by the coherent state expec- 
tation values as in subsection IV. A. Eq. (57b) is more subtle as it arises both from this 
projector as well as from the EPR correlations, Eq. (52). Because of the latter, the mean 



20 



value of the L-part of the field operator, <pL, is also described by a local wave packet even 
though nothing has been specified about L-configurations. 

The main lesson from these equations is that the simultaneous specification of the 
various has introduced some spatial coherence by coupling modes which were so far 
independent. To obtain these spatial properties, we re-use the wave packet given in Eq. 
(25). The R-wave is then given by Eq. (27) which is maximum along the classical light-like 
trajectory of Eq. (28). More interesting is the partner wave, the L component. Using Eq. 
(A6), one has 

Li Vr = -vN J cPk-^e-^ 1 ^ e -ik(x-xo) e -ifc(^H»-4^.) e -^ + c ^ _ (58) 

Two interesting properties should be discussed. First, using the stationary phase condition 
one determines the partner's trajectory defined in Eq. (34). As expected one verifies 
that the partner propagates in the opposite direction. More importantly, it is separated 
from the detected wave by the 'universal' distance given in Eq. (35). Secondly, the 
phases of the two waves are opposite when evaluated at the centers of the wave-packets, 
compare Eq. (27) and (58). This phase opposition is particularly clear when working in 1 
dimension. In this case, for any right moving wave packet, one has 

0L,v fl (x; 4>) = 4>r,v r {* - *R + x L ; 4> + vr) . (59) 

In Section IV, we have seen that a similar result obtains in 3 dimensions when working at 
the saddle-point approximation. This phase opposition originates from the coherence in 
the in- vacuum, see Eq. (14) and Eq. (A6). (Notice that it also follows from the neglect of 
the decaying mode). It has important physical consequences. It implies that the partner 
of a local Newton (Bardeen) potential dip is a local hill. For adiabatic perturbations, it 
means that the partner of a hot region is a cold region. It is to be emphasized that these 
correlations are valid for every configurations specified by Eq. (43) and not only in the 
mean. 

In brief, by having isolated from the in vacuum the (k, v) R configuration centered 
around xo at 770, we obtain a causally disconnected L configuration which is centered 
around xo — 2(770 — 2rj r ) 1^, has momentum — k, and which has the same amplitude v and 
opposite phase. Notice that these results follow from the fact that (4>}v R given in Eq. 
(56) is in fact a single wave packet of in modes. Notice finally that we have reached these 
results by making use of the complete projector Hv R which specified the amplitudes of 
all .R-modes. However these results can also be obtained when performing only a partial 
selection which leaves unspecified the amplitudes of all R-wsve packets orthogonal to the 
chosen one. The proof is given in the Appendix D. It rest on the Gaussianity of the 
distribution. 

VII. CONCLUSIONS AND DISCUSSIONS 

In the second part of this paper we have shown the following results. When considering 
(long after the reheating) the set of final configurations of a quantum field in an inflationary 
model, the detection of a coherent state describing half the modes (e.g. a right moving 
configuration) is a random process, see Eq. (49). Second, the fact that modes have been 
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amplified in pairs implies that the "reduced state" which follows from this detection is 
also a coherent state, see Eq. (52). Therefore, when specifying that some local waves have 
been detected, one introduces spatial correlations which possess definite properties. In 
particular these correlations always have the same dipolar structure since the amplification 
process is scale invariant. Moreover, the two waves in each pair are in phase opposition 
because of Eq. (A6) which tells us that the linear k dependence of the phase of Zk will 
only induce a translation of the partner's wave with respect to the chosen one. Indeed, 
this linear dependence implies that the separation between the two waves is always given 
by twice the Hubble time multiplied by the speed of the waves, in a direction specified 
by their wave vector, see Eq. (34). The phase opposition means that the spatial profile 
of the partner wave is the symmetrical of that of the chosen wave, see Eq. (59), up to a 
question of the spreads in the perpendicular directions, see Eqs. (38). It should be also 
pointed out that there also exists a strict correlation in the amplitudes of the two waves 
in each pair. This correlation in amplitude cannot be seen from the simpler treatment 
based on the two-point function since the mean has been taken before having applied the 
wave packet transformation. In addition to this discrepancy, when using a product which 
is insensitive to the sign of the wave velocity, as in Eq. (29), we obtain the three folded 
structure since the contribution of two pairs are isolated by computing the spatial overlap, 
see Figure 1. 

One might finally question if it would be possible to observe these structures. Let us 
briefly mention the different aspects which should be confronted. 

One should first analyze the statistical basis for the identification. We first notice that 
Eq. (31) results from having taken an ensemble average (which could be thought to be 
either quantum or stochastic). This mean value could be reached observationally if a 
sufficiently large number of (independent) pairs are considered. In this, one should exploit 
the isotropy and the homogeneity of the distribution, as for the temperature anisotropy 
multipoles. When dealing with modes with sufficiently high wave vectors (corresponding 
to angles smaller than a degree), this condition can probably be met. 

Second, we have access only to a portion of the configurations at recombination time, 
namely the Last Scattering Surface S, the intersection of r/ = r\ Tec with our past light 
cone, see Fig. 3. Since all pairs propagate on the particle horizon of the locus of birth, at 
recombination, they fall into three classes. First there exist pairs which do not intercept the 
LSS, such as the pair 1. These do not contribute at all to the temperature anisotropics. 
Second there exist pairs such as 2 for which only one member crosses the LSS. These 
contribute incoherently (i.e. as if = 0) to the temperature anisotropics. Third one 
finds the pairs such that both members live on the LSS. They contribute coherently to the 
anisotropies. Hence only these are responsible for the dip in the function C{9) mentioned in 
[6]. These pairs have their wave vector tangent to S. Their number is therefore limited by 
these geometrical constraints. To quantify the percentage of such pairs, one must consider 
the depth of the LSS. These aspects will be presented in a forthcoming publication. 

Moreover, the temperature anisotropies do not arise solely from the density fluctuations 
on the LSS. For a description of the various contributions, we refer to [7, 8]. Notice that the 
Doppler effect does not affect the temperature fluctuations which propagate longitudinally 
with respect to the LSS. Instead the contribution of secondary anisotropies will lower the 
level of coherence of temperature anisotropies. 

Finally, would be very interesting to have access to the velocity field on the LSS in order 
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to be able to suppress the doubling of the partners. Maybe a clever use of the polarization 
spectra might allow to reconstruct this field. Before trying to do so one might first look 
for a statistical identification of three folded structure. We re-emphasize the statistical 
character of this structure which results from an averaging procedure over pair creation 
events which individually form local dipoles. 



FIG. 3: A 2 + 1 dimensional space-time diagram in conformal coordinates (rj,x,y). The Last 
Scattering Surface S is here represented by a circle defined by the intersection of our past light- 
cone C~ with the 2d-spacelike surface z = cte,rj = rj^ss- AH pairs are created at reheating rj re h 
and propagate on light-cones. Given that the recombination is almost instantaneous, only few 
pairs, such as 3, are such that both particles intercept S. 
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APPENDIX A: THE BOGOLIUBOV TRANSFORMATION FROM INFLATION 

TO THE ADIABATIC ERA 

Since we use quantum settings, we need the Bogoliubov coefficients relating positive 
frequency modes during inflation (taken here to be for simplicity a de Sitter period) to 
modes during the radiation and matter dominated eras. 

In these three periods, the scale factor is given by 
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where rj r and r\ eq designate respectively the end of inflation and the time of equilibrium 
between radiation and matter. The shifts in the parenthesis in the second and third lines 
are necessary to parametrize the three periods by a single conformal time ??. The sudden 
pasting of these periods are such that the scale factor and the Hubble parameter are 
continuous functions. 

The positive-frequency solutions of Eq. (4) (corresponding to gravitational waves, 
the modes for density fluctuations having a different dispersion relation, see the second 
footnote) are 



^{r ] ) = ^=[l-^)e- ik \ for - oo < r, < Vr , (A2a) 



^(r,) = -L e-^ , for Vr < r, < Veq , (A2b) 



1 

V2k V Hv - Vm) 



^ at (v) = ~j= 1 - 77 v e- ikr > , for Veq < r, , (A2c) 



where we note r\ m = 4i] r — i] eq . The Bogoliubov coefficients between these modes are 

Cfa) = <~ rad <f>k ad (ri) + Pr rad * fk ad *(v) , (A3a) 
(f) r k ad (ri) = a r k ad - mat <f>T at (n) + (3 r k ad ~ mat * C a * *(v) • (A3b) 

The coefficients and (3k of the first line are given by the Wronskians 



in—rad ( A.R ain—rad* _ ( iR* ±in\ I \ A\ 

) Yk ) > Pk — ~ \Yk > Yk ) > l A4 J 



evaluated at transition time ?? r since modes satisfy different equations in each era. Similar 
expressions evaluated at the equilibrium time r\ eq hold for the coefficients between <p k ad 
and ^™ a *. One gets 

«r md = l-I-7--^L, (A5a) 



kr] r 2k 2 i] 2 



e ~2ikr) r 



2k 2 7], 



2 ' 



(A5b) 



IT 

rod-mat = j + * _ 1 ( A5c ) 

.. , 07. ... 7m )2 



k(Veq-rim) 2k 2 (r/ eq - T] r , 
-2ikrj e 

) 2h 2 (r\ eq — T]m)^ 



% ad ~ mat ) =^ ^2" ( AM ) 



These Bogoliubov coefficients have been calculated for gravitational waves in [22, 23]. 
Their results agree with ours up to constant phases which can be gauged away by a mode 
redefinition. For density fluctuations similar expressions are obtained when using the 
appropriate frequency. 

For relevant modes, i.e. modes which contribute to visible anisotropies in the CMB, 
their wave length obeys k/Ha(r] r ) = k\r] r \ <C 1. In the limiting case, to order (A;|?7 r |) 3 , one 
has 

Pk/a* k = z k = -e 4 ^ . (A6) 
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Thus the in modes during the radiation dominated era read 

m ^ rv 2,, '' +o( H' <A7> 

One then verifies that the physical modes = </>l n /ci are constant (and = H/(2k 3 ^ 2 ) 
until they re-enter the Hubble radius, i.e. when k(rj — 2r/ r ) approaches 1. This guarantees 
a scale invariant spectrum fc 3 |^| 2 = H 2 Finally, since is constant the approximation 
which consists in dropping terms of 0((kr] r ) 3 ) in Eq. (A7) corresponds to the neglect of 
the decaying mode, see the discussion after Eq. ((20)). 



APPENDIX B: COHERENT STATES FOR A REAL HARMONIC OSCILLATOR 

This appendix aims to give a self-contained presentation of coherent states, with special 
emphasis on properties which shall be used in the body of the manuscript. For more details, 
we refer to [26-28]. 

They are several equivalent ways to define coherent states. The definition we adopt 
[27] is as an eigenstate of the annihilation operator: 

a\v) = v\v) , (Bl) 

where v is a complex number. The development on this state on Fock basis is 

\v) = e-—^2^\ n )^ ( B2 ) 



n=0 

where the exponential prefactor guarantees that the state is normalized to unity {v\v) = 1. 

The first interesting property of coherent states is that they correspond to states with 
a well defined complex amplitude v. Indeed, by definition (Bl), the expectation values of 
the annihilation and creation operators are 

(v\a\v) = v , (v\a/\v) = v* . (B3) 

It is to be stressed that the variances vanish: 

Aa 2 = (v\a 2 \v) - (v\a\v) 2 = 0, Aa f2 = (v |a t2 |t>) - (v\a^\v) 2 = . (B4) 

Moreover the mean occupation number is 

(v\a^a\v) = \v\ 2 , (B5) 

in agreement with the mean value given by the Poisson distribution (B2). 

^From these properties one sees that the expectation values of the position and mo- 
mentum operators (in the Heisenberg picture, with h = 1) 

q(t) = — , p(t) = -iJ-tae-"" - tie"*) , 
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are 



q(t) = (v\q(t)\v) = -^(ve-' tujt +v*e lujt ) = J -\v\ cos(ut - <f> v ) , 
\/2lo V w 



p(t) = (v\p(t)\v) = -iJ^{ve- iwt - v*e iu)t ) = -y/to\v\ sm(ut - <j> v ) = d t q{t) . (B6) 

We have used the polar decomposition v = \v\e l< ^ v . These expectation values have a well 
defined amplitude and phase and follow a classical trajectory of the oscillator. This is 
due to the "stability" of coherent states which is better seen in the Schodinger picture. If 
the state is a coherent state \v) at a time to, one immediately gets from (B2) that at a 
later time t, the state is a coherent state given by \v(t)} = Ive'^-^) . Notice that the 
variances of the position and the momentum are 

A4»-i;, V = |. (B7) 

They minimize the Heisenberg uncertainty relations and are time-independent. Hence, in 
the phase space (q,p), a coherent state can be considered as a unit quantum cell 2nh in 
physical units (see also (C8) for the measure of integration over phase space) centered on 
the classical position and momentum of the harmonic oscillator (q(t),p(t)). In the large 
occupation number limit \v\ S> 1, coherent states can therefore be interpreted as classical 
states since Aq/q = Ap/p oc l/\v |. This is a special application of the fact that coherent 
states can in general be used to define the classical limit of a quantum theory, see [28] and 
references therein. 

One advantage of coherent states [26] is that the calculations of Green functions re- 
sembles closely to those of the corresponding classical theory (i.e. treating the fields not 
as operators but as c-numbers) provided either one uses normal ordering, or one considers 
only the dominant contribution when \v\ S> 1. In preparation for the calculations with 
field we compute the Wightman function in the coherent state \v) 

G v (t,t') = (v\q(t)q(t')\v) 

= (: mm :)v + (B8) 

where we have isolated the contribution of the vacuum. The normal ordered correlator is 



order \v\ 2 : 



(■ q(tW) :) v = -Re Ua 2 )^^ + {cf a)^'^ 

LO L 

2 

= — \v | 2 cos(w£ - <j) v ) cos(ujt' - <p v ) = q{t) q(t') . (B9) 

We see that the perfect coherence of the state, namely |(aa)„| = {a)a) v is necessary to 
combine the contributions of the diagonal and the interfering term so as to bring the 
time-dependent classical position q(t) in Eq. (B9). 

The wave-function of a coherent state in the coordinate representation is given by 

^ v (q) = Q 1/4 e -IM) 2 e iM j (Bio) 
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where v = (q + ip)/s/2uj. This follows from the definition (q\a\v) = v(q\v). ,iFrom this 
equation one notes that two coherent states are not orthogonal. The overlap between two 
coherent states is 

(v\w) = exp ^v*w — — | f | 2 — — |if| 2 ^ . (Bll) 

Nevertheless they form an (over)complete basis of the Hilbert space in that the identity 
operator in the coherent state representation {\v}} reads 




(B12) 



The measure is d 2 v = d(Rev)d(Imv). (This identity can be established by calculating the 
matrix elements of both sides of the equality in the coordinate representation {\q)}, with 
the help of (BIO).) 



APPENDIX C: TWO-MODE SQUEEZED STATES AND COHERENT STATES 

Since one can view 0k( ? ?) of Eq. (6) as the position of a complex harmonic oscillator, 
it is useful to analyze the two-mode coherent states and the two-mode squeezed states of 
a complex harmonic oscillator. Its position and momentum operators are 

1 >~ -iuit , ~t Jujt\ - ■ ■ l u 



q(t) = f + , p(t) = -^-(ane-™ 1 - d\e^) . (CI) 



Two-mode coherent states obey: 

o,r\v, R)\w, L) = v\v, R)\w, L) , &l\v, R)\w, L) = w\v, R)\w, L) . (C2) 

Therefore the expectation values of the position and momentum in the state \v, R)\w, L) 
are 

q{t) = -L(ve~ ia; * + w *e iut ), p(t) = -i\hr{ve~ iujt - w*e iujt ) . (C3) 
V2oo V 2 



As for a real oscillator, the normal ordered two-point function is given by the product of 
the mean values: 

(v, R\(w, L\ : ffltftf) : \v, R)\w, L) = q(t) q(t')* • (C4) 

A two- mode squeezed state \z)2 of this system is defined by the action of the following 
operator on the two-mode vacuum, 1 0)2 = |0, R)\0, L): 

\z) 2 = 5(r,^)|0) 2 = exp[r(e- i2 *OflO L -h.c.)] |0) 2 

= i ex p(" e+ ^ thra ^) 10)2 

= ~ l z ! 2 ) ex P l°>2 
00 

= (^lH*| 2 )£*>,fl>l«,£> (C5) 

n=0 
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where we have introduced z = — e~ l2 ^thr. The complex parameter z fully specifies the 
two-mode squeezed state. The correspondence with the Bogoliubov coefficients is made 
by z = (3 /a*, see Eq. (B13) and Appendix B in [9]. 

It is interesting to compute the projection of a two- mode squeezed state on a one- mode 
coherent state in the right sector \v, R). Using Eqs. (B2) and (C5), one gets 

f -iy s v(^) n i 



/ r>l \ 1 --\V\ ST^ V / I t\ 

(v,R\z) 2 = — e 2\ u \ ^^±-\ n: L) 




(C6) 

In the last line we have used z = (3/a* to write (1 — \z\ 2 ) as f/|a| 2 . The reduced state is 
also a coherent state. This follows from the EPR correlations in the two-mode squeezed 
state. The normalization factor can be interpreted easily. Its norm squared gives the 
probability that the system, initially in a squeezed state, is found in the coherent state in 
the right sector, irrespectively of the state in the left sector: 

Pz{v) = 2 {z\ (\v, R)(v, R\®l L )\z} 2 

\<*\ n 

\v\ 

e H 7 . (C7) 



I |2 



This gaussian probability is centered and naturally normalized to unity owing to the 
representation of unity in the coherent state basis Eq. (B12) 



We have used the decomposition v = (q + ip)/V2uj and the measure is d 2 v = 
d(Rev)d(Imv). One can see v as a stochastic variable characterized by the probabil- 
ity distribution P z {v). The variance of v is given by \a\ 2 . In the large occupation number 
limit | a | 2 S> 1 the dominant contributions of expectation values in the squeezed state \z) 2 
can be all obtained by making use of v and the distribution P z . 



APPENDIX D: THE DETECTION OF A SINGLE WAVE PACKET 

In this appendix we consider a partial projection which concern only a subset of modes. 
In fact, we shall consider the projection operator which concerns only one i?-mode and 
which acts as unity for all modes orthogonal to it. 
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1. Family of wave-packets 

Let us consider a family of positive frequency right-moving wave-packets §^(r),x.)/a(r]) 
solutions of (2). Their Fourier contend is written 



jd^Xj^^iv). (Dl) 



The functions 7^ k are parametrized by six integers, designated generically by A. Three 
of them fix the mean momentum kj = j k where the vector j = (j x ,j y ,jz) £ Z + x Z 2 
since j x > 0. The other specify the mean position x n = nx, n £ Z 3 at 770- As a concrete 
example, one can work in a box. Then the functions 7a, k are matrices. 

We assume that the & R (77, x) are orthonormal with respect to the Klein-Gordon scalar 
product 

Jd 3 x^*id v = /^7A R k7A^,k = ^aa' • (D2) 

This implies that the matrices 7^ are invertible with inverse 7^. We also assume that 
the family is complete: 

E^k7A%' = 5 3 (k-k'). (D3) 

A 

We then introduce a family of positive frequency, left-moving, wave-packets: 

/: p-ikx 
^/k^/^r- (D4) 

For reasons which shall become clear in the sequel (see Eq. (D14), (D15) and discussion 
below), we relate j L and 7^ by 

7A L ,k = -e" 2 ^7A% (D5) 

where ipk is the phase of the squeezing parameter = — \zk\e l2 ^ k . It follows that the 
matrices j L are invertible as well, and that the 3?^ are orthonormal. Since R and L wave- 
packets are orthonormal, the family {(<3?^, ^j^)} forms a complete orthonormal basis of 
the solutions of the field equation. 

Hence, the field can be decomposed as 

= E ("M + «A* A + h.c.) . (D6) 

A 

The annihilation (creation) operators are given by 

af = |^ 7 « k a«, at = |d^7A L k«k- (D7) 
and satisfy the commutation relations 

a^atf] =6 XX ,6^, (D8) 
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where fi, v stand for R, L. 

The vacuum is the tensorial product 

|0> =!]> |0, A) 2 , (D9) 

A 

where each two- mode vacuum |0, A} 2 state is defined by 

aj?|0, A) 2 = a£|0, A) 2 = 0. (D10) 
We have introduced the "tilde" tensorial product f\ to indicate that it takes into account 

A 

the indexes A which belong to Z , x Z 2 x Z 3 . 

It should be stressed that the above expression of the vacuum in terms of two-mode 
states is somehow artificial since the Bogoliubov transformation between in and out wave- 
packets will, in general, be non diagonal in A. Indeed, using Eq. (13) and (D7), one 
has: 

<R = E («AA' ~ Pxx> 8%%) • (DH) 

A' 

where 

a* xx , = jM^ai^, (3 xy = |^7jk^7A',-k, (D12) 
are not diagonal. 



2. In and out vacuum states 



To obtain the relation between in and out vacua, we start with Eq.(14), the expression 
in terms of the operators which diagonalize the Bogoliubov transformation. We then 
express the in terms of a\ so as to get 

\0in) = ^\\j^J exp (^yj J d z k z k a^ out ^ a^ out] ^j\Qout) 

= (5 j^j) ex P (E B ™ *Tn & Tt) \°° ut ) ■ ( D13 ) 

The factor L 3 /(2tt) 3 gives the density of states in a cube of size L. The matrix B\y is 
given by 

BxX ' = (S 3 /^^TaV (D14) 

In inflationary cosmology, one can fine-tune the family <I>^ so as to get a diagonal 
matrix: 

Bxx' = z x 5xx> , (D15) 
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where z\ is real and positive. One has z\ = \z^\ where k\ is the mean momentum in 
the wave packet <J>j^. Hence z\ — 1 <C 1 and n\ = nk x = 1/(1 — zf) S> 1. To obtain a 
diagonal matrix, one first fine-tunes the k dependent phase of the 7^, as specified in Eq. 
(D5). Thus one gets rid of the phase of Zk in the integrand of Eq. (D14). Second one uses 
the fact that the modulus \zk\ is a slowly varying function of k in the large occupation 
number limit when dealing with narrow wave packets in k. Then one can use Eq. (D2) 
to show that B\\r is diagonal. Let us stress that the slowly varying character of \zk\ 
follows from the large occupation number limit. Indeed, one has kSk — 1/jifc <C 1 where 
5k = dkZk — (l/rik)dk lnnfc and where we have supposed that rik is a power law. 

Using the fine-tuned L modes, the in vacuum factorizes as a product of two-mode out 
states, as in the k-basis 

l° m > = fi ( (reA+ 1 1) i /2 ex P affi 1°' A ' out h) ■ ( D16 ) 

The terms in the exponentials should be understood as a first order approximation in Sk- 
in this first order approximation, one can also treat the Bogoliubov coefficients a\\< and 
f3\\i as diagonal. In the sequel we shall work in this limit, and all equations should be 
understood as giving the leading behavior when ^ < 1. 



3. Modified correlations 

The projector on a coherent state of a out wave-packet <E>^ of amplitude v is, in total 
analogy with Eq. (42), 

U v r = \v, X,R,out)(v, X,R,out\ <g> 1 x ,l JJ <8>1a'- (D17) 

Unlike what we had with Eq. (48), this operator is unity in all sectors orthogonal to that 
defined by <E>f . 

Since the relation between in and out vacua is diagonal, we can collect the results from 
Section IV. B and V and re-express them with the A basis. The action of the projector on 
the in-vacuum is 

|i>| 2 

n„H|0m) = (n A + l) -1 e~ 2 <"A+i)| v , X,R,out) ® \(z x v*), \L,out) <g> |0, A',m)<D18) 

X'^X 

It displays coherent state correlations in that the A-th L-component is a coherent state of 
amplitude (z\v*). The probability to find the A wave-packet is given by 

1 N 

P™(„, A) = (0m|n^|0m) = — e ^+^7 . (D19) 

Notice that this amplitude is much larger than Py 1 of Eq. (48) since the coherent state 
projection concerns one mode only. (However the probability is still small because 
the resolution of coherent states is very high with respect to the occupations number 
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n x S> 1). The modified ensemble obtained by having detected § x with amplitude v is thus 
much closer to the mean than that obtained by the projector flv R . This is clearly seen by 
computing the modified expectation values Eq. (41). The 1-point functions are 



(aT R )v? = « Sx,y, ( &Zl)v? = 8 X , X , , (D20) 



or equivalently 



K!U? = v 7 i!k, ( aft,)* = 7A L k • (D2i) 



Eqs.(D21) coincide with Eqs.(53) with = wy^£ even though the projection enforced by 

TI v r is much weaker than that of Ily^ . The reason of the agreement is that the mean value 

of the modes orthogonal to <E>^ vanishes because we are in vacuum whereas it vanished 
in Eqs.(53) because we set all amplitudes to zero by the complete projection. For the 
two-point functions we have 

<4r«<^ = ^CAY™ = Km [(! - h M Wx + h,xAvf\ . (D22) 

(<r<l) v * = Km [(! - V>a/?a + h M B x \v\ 2 } . (D23) 

(D24) 

We clearly see that the 2-point functions split into two contributions. First, one finds 
the usual in vacuum expectation values in the two-mode sectors orthogonal to the chosen 
mode <&j^. Second, there is the coherent state 2-point function of amplitude v in this 
2-mode sector. 

Given these results, the conditional value of the field are 



4r,v* = = {v*T* + v **™r) , (D25a 



« = (^L,««> = {W* + z* x v nT) ■ (D25b) 



They agree with the expressions of Section V.D. because, in the mean, the expectation 
values of the unspecified modes all vanish in the vacuum. 
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